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We propose a theory of shear flow in dense granular materials. A key ingredient of the theory 
is an effective temperature that determines how the material responds to external driving forces 
such as shear stresses and vibrations. We show that, within our model, friction between grains 
produces stick-slip behavior at intermediate shear rates, even if the material is rate-strengthening 
at larger rates. In addition, externally generated acoustic vibrations alter the stick-slip amplitude, 
or suppress stick-slip altogether, depending on the pressure and shear rate. We construct a phase 
diagram that indicates the parameter regimes for which stick-slip occurs in the presence and absence 
of acoustic vibrations of a fixed amplitude and frequency. These results connect the microscopic 
physics to macroscopic dynamics, and thus produce useful information about a variety of granular 
phenomena including rupture and slip along earthquake faults, the remote triggering of instabilities, 
and the control of friction in material processing. 


I. INTRODUCTION 


Stick-slip instabilities are ubiquitous in sheared gran¬ 
ular materials ranging from pharmaceutical powders to 
earthquake fault gouge. Such instabilities are often ob¬ 
served in laboratory experiments [H-Q. At much larger 
scales, these instabilities may explain how aftershocks 
can be trigg ered by acoustic waves from earthquakes else¬ 
where |10l - ll5l| . In the context of industrial processing, 
stick-slip could result in structural damage and problems 
in mixing, causing major problems in the quality control 
of products |l6l4l8l |. However, only very few of the many 

S iin-scale numerical simulations to date, such as |1^ 
, exhibit stick-slip behavior. Most other models do not 
account for the extended elastic interaction between the 
slipping region and its surroundings. It seems that we do 
not yet know what basic ingredients of granular models 
are necessary for the apparent rate-weakening behavior 
that produces the observed unstable behaviors. For both 
practical and theoretical purposes, we need a predictive, 
first-principles description of these phenomena. 

In this paper we propose a nonequilibrium theory of 
stick-slip instabilities in dense granular flows. Our the¬ 
ory accounts for the effect of friction between colliding 
particles, and for the effects of both internally and ex¬ 
ternally driven acoustic vibrations. Our key findings are 
that interparticle friction is crucial for stick-slip insta¬ 
bilities, and that vibrations can amplify or reduce the 
stick-slip amplitude, or suppress stick-slip altogether, in 
different parameter regimes. A key ingredient of our anal¬ 
ysis is the idea that the nonequilibrium states of a driven 
granular material are characterized by its compactivity 


A = 



(l.I) 


or equivalently, its effective disorder temperature T^s = 
pX |22| - k^ . Here, V is the extensive volume of the sys¬ 


tem, Sc is the configurational entropy pertaining to slow, 
configurational degrees of freedom alone, and the A^’s 
are the internal variables, such as the density of defects 
and misalignments, that specify the configurational state 
of the granular material. The compactivity X is con¬ 
trolled by how the system is driven, and in turn controls 
the volume. In fact, X increases monotonically with the 
volume V (if the grain mass is held constant), and de¬ 
creases with increasing packing fraction. Increasing the 
shear rate would tend to dilate the system, while if one 
starts from a somewhat loose granular packing, increas¬ 
ing the vibration intensity would reduce the volume. This 
effective-temperature description is made possible by the 
fact that the slow, configurational degrees of freedom of 
the granular medium are driven out of equilibrium with 
the fast, kinetic-vibrational motions, so that the two sets 
of degrees of freedom are only weakly coupled to each 
other. 

The basic idea is that the coupling /C between the con¬ 
figurational and kinetic-vibrational subsystems consists 
of additive contributions of the fluctuations T induced 
by shearing, ^ induced by interparticle friction, and p in¬ 
duced by external vibrations (when switched on). In our 
analysis [s^ of a series of experiments by van der Elst 
et al. [^ . these are referred to as mechanical, frictional, 
and vibrational noise strengths. The effect of ^ and p on 
granular rheology is clearly illustrated there. 

In the van der Elst et al. experiments, the authors mea¬ 
sured the steady-state shear band thickness, or volume, 
as a function of the shear rate, for both spherical glass 
beads and angular sand particles in a cylindrical shear 
cell. The circular plate on top of the shear cell was ro¬ 
tated so that the shear band formed immediately under¬ 
neath the plate, and contained roughly a dozen particles 
across its width. The apparatus was connected to a trans¬ 
ducer which, when turned on, vibrated the grains at a sin¬ 
gle frequency and amplitude. Van der Elst et al. found 
that for unvibrated angular sand, the steady-state shear 
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band thickness decreases as a function of increasing shear 
rate, at intermediate shear rates, but increases with shear 
rate in the fast, inertial flow regime. This is in contrast 
to smooth spherical glass beads for which the shear zone 
thickness increases monotonically with the shear rate. 
Moreover they found that external acoustic vibrations 
cause compaction in both angular sand and smooth glass 
beads at slow shear rates. Their experimental measure¬ 
ments and our quantitative fit to their data are shown 
in Fig. 1 in our earlier analysis [s^. In Fig. [T] below, 
we also plot the dimensionless compactivity, instead of 
the volume or shear band thickness, as a function of the 
shear rate, albeit at a higher confining pressure than in 
the van der Elst et al. experiments. 

Our analysis has shown that acoustic vibrations flu¬ 
idize the granular material; that is, vibrations unjam the 
granular material and significantly reduce the flow stress 
at slow shear rates. In addition, frictional dissipation 
between particles suppresses shear dilatancy at interme¬ 
diate shear rates, and is responsible for the nonmono¬ 
tonic variation of the steady-state shear band thickness 
with shear rate seen in those experiments. This is shown 
by a direction comparison between the dashed and solid 
curves for frictionless and frictional particles in Fig. [T] 
Nonmonotonic rheology is suggestive of inherent insta¬ 
bilities [13, stick-slip failure being one example. This 
motivates our present theoretical analysis. We show in 
this paper that interparticle friction is indeed responsible 
for stick-slip at intermediate shear rates and high enough 
pressure. In addition, vibrations can amplify or reduce 
the stick-slip amplitude, or suppress stick-slip in different 
parameter regimes. 



FIG. 1: (Color online) A key ingredient for the emergence of 
flow instabilities is the nonmonotonic variation of the steady- 
state compactivity as a function of the imposed shear 
rate 7 . (The quantity y®® increases monotonically with the 
layer thickness or volume, so the layer thickness follows the 
same qualitative behavior.) For clarity, this is shown here for 
a pressure p = 10 MPa, which lies below the critical pressure 
necessary for stick-slip to occur (see Fig. [S] below). It is im¬ 
plicitly assumed that physical insight gained from laboratory 
experiments at lower pressures can be extrapolated to higher 
pressures in the MPa range. 


The rest of the paper is structured as follows. Sec¬ 
tion El is an exposition of the thermodynamic theory of 
driven granular media. There, we demonstrate how the 
first law of thermodynamics allows us to infer the man¬ 
ner in which the compactivity X evolves as a function 
of time in sheared and/or vibrated granular media. This 
is followed by Sec. IIIII in which we decouple the effect 
of different driving mechanisms, and specify the config¬ 
urational state of the granular media when either shear¬ 
ing and vibrations but not both is present. Then, in 
Sec.HVl we discuss the microscopic model of shear trans¬ 
formation zones (STZ’s) [H, [H, [sJl which provide 

a direct connection between mechanical stress and irre¬ 
versible, nonafhne granular rearrangement that accounts 
for dense granular flow. Upon completion of the theo¬ 
retical development, we show in Sec. |V] the emergence of 
stick-slip from our model, examine the effects of external 
vibrations at fixed frequency and amplitude, and discuss 
the conditions under which stick-slip instabilities are am¬ 
plified or suppressed. We conclude our paper in Sec. EH 
with a discussion of open questions and future directions. 


II. NONEQUILIBRIUM THERMODYNAMICS 
OF DRIVEN GRANULAR MEDIA 

The following theoretical development is almost the 
same as, but not identical to, that presented in Ref. [^ . 
Here, we explicitly account for the work done by exter¬ 
nal acoustic vibrations - which need not be negligible 
- and clarify the nature of the thermal temperature T. 

We also discuss how external vibrations directly control 
the configurational degrees of freedom of the system in 

Sec. EH 

Consider a granular medium in contact with a heat 
bath at temperature T Ki 0 (we set T = 0 later but 
retain the variable T at present), with corresponding en¬ 
tropy St- Then, when the system is driven by a shear 
stress s in the presence of a pressure p, and if vibra¬ 
tions do work on the medium at a rate pVY/r (where 
the time scale r is conveniently chosen to be the inertial 
time scale, discussed in Sec. 113 below, and we have fac¬ 
tored out pV to make Y dimensionless), the first law of 
thermodynamics reads 

TSt = 

r 

= VsY^-pXSc-p'^(-^'] Aa-b^^^(2.1) 

where 7 P* is the plastic shear rate. As in Eq. dm, the 
Ac’s are the internal variables, such as the density of 
flow defects, that specify the configuration of the granu¬ 
lar medium. Note that while T ft; 0, the quantity TSt 
need not be small, as energy is being dissipated into the 
environment. On the other hand, the second law of ther¬ 
modynamics says that the rate of change of total entropy, 
being the sum of the thermal and configurational contri- 
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buttons, must be nonnegative: 

S = Sc + St>0. (2.2) 

Rearranging Eq. m to obtain an expression for Sc, 
substituting into Eq. (IQ) . and using the fact that each 
individually variable term in the resulting inequality 
must be nonnegative , we arrive at the second-law 

constraints 


W = Vs'jP^-pY^ 


f dV 


\dKo 


vVY 

Ac -b - - > 0; (2.3) 


Sc 


(pX - T)St > 0. (2.4) 


In arriving at these two constraints, we have arranged 
terms in such a way that terms pertaining to the de¬ 
grees of freedom that belong to the same subsystem are 
grouped together. The dissipation rate W, as defined 
in [36l - l38l| , is the difference between the rate at which in¬ 
elastic work is done on the configurational subsystem and 
the rate at which energy is stored in the internal degrees 
of freedom. The second constraint implies that pX — T 
and St must carry the same sign if they are nonzero, so 
that 


TSt = -VIC{X, T){T-pX)= Q, (2.5) 

where /C(X, T) is a non-negative thermal transport coeffi¬ 
cient. It is already clear from this analysis that p X plays 
the role of a temperature. A heat flux Q flows between 
the granular subsystem and the reservoir when the two 
subsystems are not in thermodynamic equilibrium with 
each other. 

We can now set T = 0, and define the dimensionless 
variables, to be used in the rest of the development, as 
follows: 




( 2 . 6 ) 


Thus X is the dimensionless compactivity, and p is the 
shear-stress-to-pressure ratio, vz is the excess volume of 
STZ’s - loose spots where irreversible particle rearrange¬ 
ments occur, incurring local topological change. The ex¬ 
tensive volume V of the granular medium is a sum of the 
total volume of grains Vb plus the contributions of the 
configurational degrees of freedom, such as STZ’s and 
misalignments. In the simplest approximation, it varies 
linearly with y: 


and replacing Sc by x through the chain rule of differ¬ 
ential calculus, the equation of motion for x becomes 


eiX = Mi''"' + — - ICix)x- 

T 


( 2 . 8 ) 


The next step is to determine the transport coefficient 
JC{x) in Eq. (12.81) . Because /C(x) couples the configu¬ 
rational and kinetic-vibrational subsystems, we assume 
that it consists of additive contributions from shearing, 
interparticle friction, and vibrations. Our strategy is 
to first determine JC{x) for frictionless grains, then add 
in the effect of friction afterwards. We know that in 
the absence of vibrations and interparticle friction, the 
steady-state compactivity is uniquely determined by the 
dimensionless shear rate q = that is, x'*’* = x('Z)- 
As such, direct substitution into Eq. (1^ in this case 
yields fC = W/{Tx{q)). Here the dimensionless work 
rate W = r/iyP* is proportional to the mechanical noise 
strength T, to be discussed in greater detail in Sec. ITVl 
below. 

On the other hand, if the granular medium is vibrated 
but not sheared, the steady-state compactivity is de¬ 
termined by the vibration intensity or noise strength 
p m, |i3|: x'''' = X(p)- As such /C = Y/{tx{p))- The 
quantity p should be a function of the vibrational fre¬ 
quency and amplitude; past studies [3, [3) IHl suggest 
that p is proportional to the vibrational amplitude, and 
the square of the frequency. However, the detailed de¬ 
pendence is unimportant for the purposes of this paper 
which, as in the van der Elst et al. experiments 3l|, fo¬ 
cuses on one single amplitude and frequency, unless oth¬ 
erwise specified. 

If the granular material is sheared and vibrated at the 
same time, it is plausible that the shear rate sets a time 
scale below which vibrations cannot compete with shear¬ 
ing in allowing the system to explore different configura¬ 
tional states by rearranging constituent grains. Referring 
to Fig. [U and foreshadowing the discussion in the next 
subsection, the effect of vibrations to increase the pack¬ 
ing fraction relative to an unvibrated granular medium 
diminishes above a certain shear rate. Because the in¬ 
tensity of vibrations is proportional to the square of the 
frequency [3, US] : we propose that the diminution of the 
effect of vibrations is determined by the ratio q'^/p. Thus, 
for smooth frictionless grains, the transport coefficient is 


/C(x) = - 

T 


w 


Y 


.X(9) X(p). 


(2.9) 


X = l + ei(x-Xr), (2.7) 

Xo 

where ei can be interpreted as an effective volumetric 
expansion coefficient, and Xr denotes the compactivity 
of a reference state. 

The equation of motion for x is derived from the ex¬ 
pression for Sc- To this end, eliminating St in Eq. (HID 
using Eq. (1^ for an expression for Sc, dividing by pVg, 


where r = exp(—g^/p) is the relative weight of the con¬ 
tributions from vibrations versus shear. Note that this is 
for frictionless grains. Of course, r — Q when the granu¬ 
lar medium is not vibrated, i.e. p = 0. In the cases when 
the granular medium is subjected to either shear or vi¬ 
brations but not both, Eq. (12.91) reduces to the results 
discussed in the preceding paragraphs. In analogy to the 
direct proportionality between the mechanical noise T 
and the plastic work of deformation W alluded to above. 
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we propose that the noise strength p associated with vi¬ 
brations is proportional to the work of external vibrations 
itself: Y = Aqp for some constant Aq. 

For the case of angular, frictional grains, we propose 
to modify the coupling coefficient K. as follows: 


nx) 


1 

T 


'W + F Y ' 

. xiq) ^^xip)\' 


( 2 . 10 ) 


where F denotes the dissipative effect of friction. As 
such, the equation of motion for y, Eq. (15.31) . becomes 




w + r 

T 


1 

T 


'W + F Y ' 

. x(q) ~^^x(p). 


X- 


( 2 . 11 ) 


The way F appears in the equation ensures that at large 
shear rates, the dilatational effect of shearing trumps the 
compaction brought about by external acoustic vibra¬ 
tions and frictional interaction, as suggested by experi¬ 
ments such as [3l|. This can be seen directly by setting 
X = 0 in Eq. (j2.11l) and solving for the steady-state com- 
pactivity y®®. In analogy to W and E, we propose that 
F is proportional to what we term the “frictional noise 
strength” to be discussed in Sec. |Vl 


III. SHEAR-RATE DILATION AND 
VIBRATION-DRIVEN COMPACTION 


At this point let us discuss the quantity x(g), which de¬ 
notes the disorder temperature of a driven granular mate¬ 
rial in the absence of vibrations and interparticle friction. 
One example of a frictionless granular material is an ide¬ 
alized hard-sphere system. We pointed out in that 
this quantity approaches some constant yo in the limit 
of small shear rate q, and becomes a rapidly increasing 
function of q once the shear rate becomes comparable to 
the rate of intrinsic structural relaxtion, thereby implying 
shear-rate dilation. A convenient way to interpolate be¬ 
tween these two limiting regimes is inspired by the Vo gel- 
Fulcher-Tamann (VFT) form in glass theory |4^.TS|| 
for the inverse function q(x)- 


- = — exp — + Q!eff(x) , 

q qo lx 


where 


aeff(x) 


Xi 


X-Xo 


exp 


XA-XeJ 


(3.1) 


(3.2) 


In the rest of the paper, we have chosen yo = 0.3, 
yi = 0.02, XA = 0.33, A = 2 and go = 2; these values 
are identical to those used in to analyze experimen¬ 
tal data for auto-acoustic compaction |3lj | at intermedi¬ 
ate shear rates. (The appearance of the constant 3 in 
Eq. (13.21) is purely a matter of convention.) 

The choice A = 2 in Eq. (ED stipulates a rate¬ 
strengthening response of the granular material to ap¬ 
plied shear, at least in the large q regime. An earlier 


analysis by Daub and one of us [IJ involved a model 
for which A < 1, which corresponds to rate-weakening 
behavior; stick-slip instabilities did occur in that case. 
However, experience shows that hard-sphere systems, or 
“barebone” granular materials, should be rate strength¬ 
ening [ 2 ^. A key result of this paper is that the con¬ 
dition A < 1 can be lifted in the presence of other mi¬ 
croscopic mechanisms, such as interparticle friction, that 
alter the hard-sphere rheology, allowing for the possibil¬ 
ity of stick-slip. In other words, we isolate the mechanism 
for rate-weakening behavior and stick-slip instabilities in 
a naturally-occurring granular medium. 

Meanwhile, the quantity y(p), which denotes the 
steady-state compactivity of a vibrated granular 
medium, should be a decreasing function of increasing 
vibration intensity p. (As we pointed out above, p is 
a function of the vibration frequency and amplitude.) 
This is because an increase in vibration intensity causes 
the grains to explore more pos sible configurations and re¬ 
duce the packing fraction , dissipating energy in 

the process. Since we consider only one vibration inten¬ 
sity, we do not propose a functional form for y(p), but 
stipulate that it takes on one of two values depending 
on whether the fixed-intensity vibration is turned on or 
off. The combined effect of shearing, vibrations, and fric¬ 
tion on the steady-state compactivity y®® is illustrated in 
Fig. m above for vibrated and unvibrated granular media 
composed of frictional and frictionless particles. 


IV. MICROSCOPIC MODEL OF STZ’S 

We now turn to the Shear-Transformation-Zone (STZ) 
theory of dense granular flow [H, [H, HI, which 
directly connects macroscopic dynamics to the grain- 
scale physics of deformation. The latter refers to the 
irreversible, nonafhne granular rearrangements, along 
with noise induced by friction and external vibration. 
The STZ theory has been invoked in our successful 
effort to reconstruct and explain constitutive friction 
laws [^.Idslld^. 

Recall our physical picture that in dense granular flow, 
irreversible particle rearrangements occur at rare, non¬ 
interacting soft spots with excess free volume known as 
STZ’s. The applied shear stress defines a direction rel¬ 
ative to which STZ’s can be classified according to ori¬ 
entation, with total numbers N+ and respectively. 

Upon application of shear stress in the “plus” direction, 
STZ’s of the minus type easily deform to become plus- 
type STZ’s. However, plus-type STZ’s rarely flip and ac¬ 
quire the minus orientation; rather, they are annihilated 
readily by noise. This is described by a master equation 
of the form 

tN± =n{±p,x)N^-n{TP.x)N±+f . 

(4.1) 

Here, we assume that the attempt frequency F = F -|- 
p, being the sum of mechanical and vibrational noise 
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strengths, does not include contributions from be¬ 
cause interparticle friction serves to dissipate energy but 
does not open up or close voids. The inertial time scale 
T = a\/ Pg/p 143) with a, pc, and p being the char¬ 
acteristic grain size, grain material density, and confin¬ 
ing pressure, is the typical duration for a pressure-driven 
granular rearrangement. 

If N is the total number of grains, the plastic strain 
rate is 


and (lO) . do not contain the small factor A, but Eq. (14.81) 
does. This provides the basis for a further simplification 
- the quasisationary approximation - under which we can 
replace A and m in the equations by their quasistationary 
values A®*^ and rrf'^. We have argued that A®*^ = 

To calculate invoke the Pechenik hypothesis, which 
states that the mechanical noise strength T is propor¬ 
tional to the plastic work of deformation dissipated per 
STZ gil: 


^ x)N- - 7^(-/i, x)N+] . 


Introducing the intensive variables 


A = 


N+ + N_ _ 
N ’ 


m = 


N+-N_ 

N+ + N_' 


(4.2) 


(4.3) 


which denote the density and orientational bias of STZ’s, 
and analogously the rate factors 


C(m,x) = ^ (^(M)X)a)) ; 

T(n vi = ^(m,x) -^(-A^)X) 

n{p,x)+n-^^,xy 


(4.4) 

(4.5) 


Eqs. (14.11) and (14.21) can be rewritten in terms of intensive 
quantities alone: 


tA = f(A®'i-A); (4.6) 

rm = 2C{p,x)iTip,x)-m)-fm-T^m] (4.7) 
= 2eoAC{p,x)iTip,x) -m), (4.8) 


Note that A and m are not the only internal state 
variables; others include the orientational bias of angular 
grains [^. However, at least for grains with aspect ratio 
close to unity, those variables are not coupled to Eq. (14.81) , 
and therefore play no role in plastic shear deformation. 
Thus, internal variables other than A and m are expected 
to be immaterial in the present analysis, in which we are 
not explicitly concerned with the volumetric response of 
the system. 

We now substitute Eqs. (g^, dUZl), and (Oil into 
Eq. (12.311 for the dissipation rate W which, according 
to the second law of thermodynamics, must be non¬ 
negative. We then arrive at several independent con¬ 
straints on the dynamics of the system, according to 
the procedure outlined in [H, 113, l35l - l38l| . One of these 
constraints says that the STZ density at equilibrium 
equals A®‘i = 2e“^/^. Another tells us that T{p,x) = 
tanh(eoAt/ezx)) where ez is the excess volume vz of each 
STZ normalized by the grain size a^. 


A. Mechanical noise and steady-state dynamics 

Since STZ’s are rare and non-interacting, A <§; 1. Ob¬ 
serve that the equations of motion for A and m, Eqs. (gH) 


eoAtoA®q 


‘^Rq(T{p,x) - my 


(4.9) 


In the last equality we replaced C{p, x) by some constant 
i? 0 ) for the STZ flipping rate should not depend strongly 
on the shear stress as long as the shear-stress-to-pressure 
ratio p is less than unity. Then the stationary version of 
Eq. (I4.7|l reads 


2i?, 


(j'{p,x) 


— m 


I - 


mp 

Po 


— mp = 0, (4.10) 


and the stationary value of m is given by 




2p 


1 + y-T{p,x )+^ 

Po 


Mo 

'2p^ 


1 + -^T{p,x)+ -jy- 

Mo ^-^0, 


1 2 


-^Tip,yfi.n) 

Mo 


In particular, when p = 0, we have 

r^eq ^ \T{p,x), if (m/mo)7'(m,x) < 1; 
\mo/m, if (m/Mo) T(m, x) > 1- 


(4.12) 


Thus, in effect, po sets the minimum flow stress of the 
system in the absence of external vibrations. However, 
when p ^ 0, the system is un-jammed, and flows at arbi¬ 
trarily small shear stress s. 


B. Yield stress parameter and apparent weakening 

The yield stress parameter po, which originates as a 
proportionality parameter between the mechanical noise 
strength T and the plastic work of deformation, plays a 
crucial role in the present analysis. In we assumed 
that Po is a rapidly decreasing function of the kinetic- 
vibrational temperature below the glass transition tem¬ 
perature. Yet it could equally well be a decreasing func¬ 
tion of the compactivity x if we extend the analogy be¬ 
tween the thermal and effective temperatures. There 
are two ways to understand this. First, intuition sug¬ 
gests that the amplitude of the fluctuations generated 
by shearing decreases as the granular medium becomes 
more tightly packed. Thus po should decrease with in¬ 
creasing y. Second, one expects that a larger shear stress 
is needed to create voids in a more densely-packed gran¬ 
ular medium in order for it to flow. In this connection. 
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we pointed out above that /zq controls the minimum flow 
stress in the absence of vibrations. Indeed, a /ip that de¬ 
creases rapidly with y when y is below the steady-state 
compactivity xo in the slow-shear limit, such as shown 
in Fig. [2 is an essential ingredient for the emergence of 
stick-slip instabilities here. 



Compactivity x 


FIG. 2: (Color online) The parameter fiQ, which provides a 
measure of the apparent yield stress, as a function of the 
compactivity x- This quantity must be a rapidly decreasing 
function of x whenever x < XOi &nd cease to be so above XOi 
as indicated by the vertical dashed line, (xq is the steady- 
state compactivity in the slow shear limit.) This condition 
conforms with the intuition that the stress needed to unjam 
a granular medium increases with the packing fraction. This, 
combined with the rapid increase of the steady-state com¬ 
pactivity x'*'* below Xo for frictional grains shown in Fig. [T] 
provides the apparent weakening that permits stick-slip in¬ 
stabilities at high enough pressure. Here we have used the 
interpolation Eq. (14.1311 with pi = 0.3 and p 2 = 3.7 x 10“®. 


On the other hand, po must cease to be rapidly de¬ 
creasing when X > XOj and be bounded from below by 
a nonzero value, so that the energy dissipated by shear- 
induced fluctuations (to open up voids and create STZ’s, 
for example) cannot exceed the energy fed into the gran¬ 
ular medium itself. Thus we use 


V. EMERGENCE OF STICK-SLIP 
INSTABILITIES 


A. Equations of motion 


The complete constitutive description consists of two 
dynamical equations, one for the shear-stress-to-pressure 
ratio p and one for the compactivity x- The equation for 
p is a statement of linearity between increments in stress 
and elastic strain. Thus, 

A = (G/p)(7-7p'), (5.1) 


where G is the aggregate shear modulus, and 7 is the im¬ 
posed shear rate. The plastic strain rate is given directly 
by Eq. ()4.8p : 

yPi = [Tifi, x) - to(p, x)] • (5-2) 

T 

Meanwhile, the equation for the disorder temperature x 
describes how the fluctuations T induced by shearing, 
and p induced by external vibrations, drive the system 
towards different configurational steady states, described 
respectively by xio) ^-nd x(p). Additionally it describes 
how friction-induced noise ^ when the material is sheared 
results in further dissipation and compaction: 


X 


2eopoe 




r 1- 




X 


X 


x{q)j x{q) 


^ Mp 

rei 


1 — exp 



X 

x(p). 


(5.3) 


To derive Eq. (15.31) from Eq. (12.111) directly, it suffices 
to identify in Eq. (j2.11|) W = TeopoA, V = Aop, and 
F = ^coMoA. The first two of these equalities have been 
discussed in connection with Eq. (I 2 J 1 ) : they imply that 
the work rates of shearing and vibrations are proportional 
to the respective noise strengths. The last equality con¬ 
nects the frictional dissipation F with the noise strength 
^ in the same manner. It simply says that frictional dissi¬ 
pation increases with the number of STZ’s present in the 
granular medium. We discuss in the following subsection 
the noise strength ^ in greater detail. 


Po = (4.13) 

to interpolate between the large- and small-x regimes. 
We argue below that a po with these qualitative fea¬ 
tures is one of two ingredients that provide apparent 
rate-weakening at intermediate shear rates. The other 
ingredient is the nonmonotonic variation of the steady- 
state compactivity x’’’’ with shear rate, shown in Eig. [T] 
Referring to Eq. (12.111) and the subsequent discussion, 
this is provided by appropriate noise strengths p and ^ 
associated with vibrations and interparticle friction. The 
apparent weakening shows up in the phase plot of shear 
stress s versus plastic strain rate 7 P* in Eig. 2] 


B. Noise strength ^ associated with interparticle 
friction 

The steady-state solution to Eq. (|5.3p shows that if 
^ varies as ( 7 ^*)^ at small shear rates, and saturates to 
some constant ^0 at large 7 ^', then it is possible for x to 
decrease below the zero-shear-rate limit xo for a range 
of intermediate shear rates. This is illustrated in Fig. [T] 
We argued in that this nonmonotonic variation of x 
explains the occurrence of auto-acoustic compaction in 
that regime. Thus we employ the interpolation 

i = tanh[(r/7P')^] 


(5.4) 
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between the two limiting regimes with rj being a time 
scale at which the compactional effect of interparticle fric¬ 
tion becomes prominent. 

At somewhat larger shear rates, x increases as a func¬ 
tion of the shear rate but remains smaller than xo- ix 
only exceeds xo and diverges at very fast shear rates.) 
Here, /tq is a rapidly decreasing function of increasing x- 
As a result, /tq decreases markedly as increases. This 
apparent rate-weakening behavior, summarized in Figs.[T] 
and [21 permits stick-slip instabilities in unvibrated fric¬ 
tional granular materials for this range of intermediate 
shear rates. One might ask whether stick-slip could also 
occur for unvibrated frictionless particles, which corre¬ 
spond to the dashed red curve in Fig. [T] for which there 
is also a rapid increase of x as a function of yP' even 
when X < Xo- However, vibrations fluidize the granular 
medium and cause weakening. That is, a vibrated gran¬ 
ular material can experience irreversible deformation at 
arbitrarily small stresses, and the behavior of /xg becomes 
immaterial in that regime. There is no physical mecha¬ 
nism that brings the steady-state x below xo and then 
up again. Therefore stick-slip instabilities do not occur 
in that case. 


C. Results 

We now present the key predictions and implications 
of the STZ granular flow model. Equations and 
(15.311 are numerically integrated using an adaptive, im¬ 
plicit time-stepping scheme, for a range of values of the 
pressure p and imposed shear rate 7 , and initial values 
fi{t = 0) = 10“^ and x{t = 0) — 0.18. We focus on 
angular quartz sand with a typical particle diameter of 
a = 350 /xm, mass density po = 1600 kg m“^, and ag¬ 
gregate shear modulus G = 109 MPa. Because of the 
complexities in the functional dependencies of the vari¬ 
ous parameters, we have chosen not to use dimensionless 
parameters in our results. Instead, we express the shear 
rate in units of s“ ^, and pressure in units of MPa, to make 
the pressure-dependent behavior more transparent. We 
take Ro = 1, Aq = 0.01, eg = 1.5, ez = 0.5, and ei = 0.3 
in Eqs. (15.21) and (15.31) . and ^g = 4 x 10“^, r/ = 0.013 s, 
and Xo = 0.3. For comparison we also compute the re¬ 
sponse of frictionless spherical beads to shear, with ^g = 0 
and otherwise identical material parameters. Referring 
to Sec. m above, our choice of x(<?) automatically stipu¬ 
lates rate-strengthening rheology in the limit of large 7 ^^ 
Unless otherwise specified, we take x(p = 0) = 0.2 in the 
absence of external acoustic vibrations, and x(p) = 0.18 
and p = 5 X 10“"* in its presence. Most of these param¬ 
eters are identical to those used in to model auto¬ 
acoustic compaction at intermediate shear rates. The 
parameter values are either determined empirically, such 
as for ei and ez, or are chosen to quantitatively describe 
the experimental findings of [^ . In addition, we choose 
Pi = 0.3 and p2 = 3.7 x 10“^, so that poix) decreases 
rapidly as a function of x when x < Xo = 0.3, and levels 


off above xo- 



0.0 0.5 1.0 1.5 2.0 

Shear strain y 


FIG. 3: (Color online) Variation of (a) shear stress s and 
(b) compactivity x with shear strain, in the presence (solid 
curves) and absence (dashed curves) of interparticle friction, 
and in the absence of external vibration. Observe that stick- 
slip occurs only in the presence of interparticle friction, which 
provides a means for the compactivity to fall below xo = 
0.3, shown with the dotted line in (b), in the ideal hard- 
sphere model, and facilitates the emergence of instabilities. 
The slip phase of stick-slip coincides with a sharp increase in 
the compactivity, or volume dilatation. 

Figure in makes transparent the essential role of inter¬ 
particle friction in the emergence of stick-slip. No stick- 
slip occurs in the absence of friction, that is, when ^ = 0 
in Eq. (15.31) . (For clarity we plot the shear stress s instead 
of the shear stress to pressure ratio p.) During the initial 
loading, the red and blue curves for stress versus strain in 
Fig.jSKa) describe almost purely elastic deformation; thus 
X remains constant at its initial value. This is indepen¬ 
dent of the presence of interparticle friction. Then, when 
the stress exceeds the yield stress for that value of x, the 
plastic strain rate 7 ^' jumps up, as does x according to 
Eq. (15.31) . A dramatic stress drop follows immediately, 
as a consequence of Eq. (EH), causing the arrest of slip. 
In the absence of interparticle friction, the compactivity 
X remains large, so that po{x) stays small. As such, the 
material remains soft, and no further slipping occurs. 

However, this may not be the case if there is fric¬ 
tion. As we argued above and in [ 3 ^, the noise ^ as¬ 
sociated with interparticle friction dissipates energy and 
cools down the system in the volumetric sense; it causes 
the compactivity x to fall below xo = 0.3, seen also in the 
solid curves in Figs. [T] and Hb). This is a regime where. 

















according to Fig. [2l fJ,o{x) is a rapidly decreasing func¬ 
tion of X, providing the necessary weakening effect for 
stick-slip instabilities to occur. Note the important role 
played by the y-dependence of the apparent yield stress 
^o(x)- Interparticle friction increases the shear stress at 
all times, conforming with the notion that friction in¬ 
creases the difficulty for slip to occur. Note also that the 
slip phase corresponds to a dramatic increase in the com- 
pactivity X; which increases monotonically with the layer 
thickness or volume [s^. A direct connection to labora¬ 
tory measurements of volume or layer thickness can be 
made once the quantitative relation between y and the 
volume is established. Thus the disorder temperature y 
provides a natural link between stick-slip dynamics, vol¬ 
umetric effects, and energy dissipation. 

Another way to visualize the dynamics is to construct 
a phase plot for the evolution of shear stress and com- 
pactivity, as functions of the plastic strain rate in 
each stick-slip cycle. This is shown in Fig. |4l where the 
imposed strain rate 7 is indicated by the vertical dotted 
line. In each case the material starts from nearly zero 
shear stress and is loaded elastically; thus yP' = 0 and 
is not shown on the phase plot. Then, when the shear 
stress exceeds the yield stress, plastic deformation sets 
in; yP' increases dramatically, depicted by the horizontal 
portions of the phase plots shooting towards the right. 
The release of plastic strain results in a dramatic stress 
drop, after which the the phase curve either converges 
to a fixed point in the case of frictionless grains, or a 
limit cycle in the case of stick-slip in frictional grains, for 
the values of pressure p and imposed strain rate y shown 
here. Creep occurs during the stick phase; that is, plastic 
deformation occurs at a much slower, but nonzero, rate 
than the imposed strain rate. The shear stress increases, 
and the compactivity decreases, so that the granular layer 
contracts, until the applied load overcomes the frictional 
resistance. This is followed by a dramatic increase of 
the plastic strain rate, or slip rate, which corresponds to 
the compactivity surge and stress drop. Because of the 
elasticity of the granular medium itself, the applied load 
“sticks” again, and the stick-slip cycle repeats itself. The 
apparent weakening behavior is evident from Fig. |31[a); 
the stick-slip phase curves always traverse the limit cycles 
in the clockwise sense. Stress drops occur at large plastic 
strain rates, while stress increases develop at lower yP*. 

Of interest is the rounding of the stress-strain curves 
near each slip event; that is, in each stick-slip cycle, the 
rate of increase of the shear stress s during the stick phase 
decreases with increasing shear strain. This creep is also 
seen in experiments [l| . Referring to Eq. (ED, the plastic 
strain rate yP' must therefore be nonzero even during the 
stick phase. This is reminiscent of preseismic slip in the 
context of earthquakes; that is, small but nonzero plas¬ 
tic deformation prior to the large slip events that corre¬ 
spond to dramatic stress drop. Evidence of this is shown 
in Fig. [51 Panel (a) shows once again the stress-strain 
curves. In panel (b), we plot the plastic strain rate yP^ 
versus the total shear strain y, and show that the plastic 



Plastic strain rate yP'(s ') 

FIG. 4: (Color online) Variation of (a) shear stress s and (b) 
compactivity y with the plastic strain rate yP* in the pres¬ 
ence (solid curves) and absence (dashed curves) of interpar¬ 
ticle friction, and in the absence of external vibration. The 
pressure values are p = 40 MPa (blue) and 60 MPa (red), 
and the imposed strain rate is y = 10^ s“^, indicated by the 
vertical dotted line. Stick-slip occurs for frictional grains, and 
each of the corresponding phase curves converges to a limit 
cycle. For frictionless grains, each of the phase curves con¬ 
verges to a fixed point, indicated by each of the filled circles, 
on the phase diagram. (The fixed points of y versus yP* for 
p = 40 and 60 MPa overlap with one another.) In the case of 
frictional grains, the applied load “sticks” left of the vertical 
dotted line, and creeps at a slower but nonzero rate than the 
imposed strain rate, while the compactivity decreases. Once 
the shear stress exceeds that needed to overcome frictional 
resistance and initiate slip, both the compactivity and the 
plastic strain rate increase dramatically, and the stress drops. 
Then the applied load sticks, and the shear stress builds up 
again during the new stick phase. Thus the stick-slip cycle 
repeats itself. This representation of stick-slip dynamics can 
be directly compared to that in Fig. [51 

strain rate is nonzero and increasing even during the stick 
phase of each stick-slip cycle. During the slip phase, yP^ 
increases very rapidly, and peaks at a value that is a few 
orders of magnitude above the imposed strain rate y. As 
a result of the stress drop, the plastic strain rate plum¬ 
mets almost instantaneously and the load sticks again. 
We plot in panel (c) the accumulated plastic strain 



as a function of total shear strain y, and show that while 
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FIG. 5: (Color online) Variation of (a) shear stress s, (b) plas¬ 
tic strain rate and (c) accumulated plastic strain with to¬ 
tal shear strain, in sheared, unvibrated angular sand, a gran¬ 
ular material composed of frictional grains. Observe from (b) 
that 7 *’* is slowly increasing during the stick phase of each 
stick-slip cycle; this accounts for the the apparent softening 
seen towards the end of each stick-slip cycle in (a). Panel (c) 
shows that while much of the plastic strain is accumulated 
during the slip phase, there is plastic strain during the stick 
phase as well. 


much of the plastic strain is accumulated during the slip 
phase of each stick-slip cycle, there is indeed a slow build¬ 
up of plastic strain, or creep, during the stick phase. This 
corroborates with the result shown in Fig. [SKb). While 
not shown here, the same qualitative behavior is observed 
in vibrated frictional sand particles when stick-slip oc¬ 
curs. 

Figures | 6 (a)| shows sample stress-strain response for 
angular sand, with external vibrations turned on and off 
alternately. Here, the imposed shear rate is 7 = 2.4 x 10^ 
s“^, and the pressure values arep = 45 and 50 MPa. Sev¬ 
eral interesting features are observed. First, at the lower 


Vibrations suppress or amplify stick—slip, y = 2.4 X 10^ s * 

60i . 


(a) 


— p = 45MPa 

— p = 50 MPa 






10 


Shear strain y 


45 


40 


35 


30 


25 


20 L 


Vibrations suppress or amplify stick-slip, y = 2.4 X 10^ s ' 

(b) ... 





-p = 45 MPa, unvibrated 

- p = 45 MPa, vibrated 
-p = 50 MPa, unvibrated 

- p = 50 MPa, vibrated 


10 


100 1000 10'* 

Plastic strain rate 7 '’’ (s “*) 


10 " 


FIG. 6 : (Color online) Vibrations suppress or amplify stick- 
slip in different parameter regimes, (a) Shear stress s a func¬ 
tion of shear displacement 7 for sheared angular sand, for 
different magnitudes of the normal stress p, at applied shear 
rate 7 = 2.4 x 10^ s“^. External acoustic vibration is turned 
on (with intensity p = 5 x 10““^) for 6 < 7 < 8 , and off oth¬ 
erwise. At the lower pressure of p = 45 MPa (green curve), 
vibrations amplify stick-slip amplitude before atttenuating it. 
At the higher pressure of p = 50 MPa (red curve), vibrations 
amplify the stick-slip amplitude and period. The red curve 
for p = 50 MPa has been offset by 10 MPa for clarity, (b) 
Variation of shear stress s with plastic shear rate 7 *"* under 
the conditions in (a) above, for shear strains 4 < 7 < 8 . The 
dashed lines denote the behavior under vibrations of fixed 
intensity, while the lighter shades indicate transients before 
approaching a limit cycle or fixed point. At the lower pres¬ 
sure of p = 45 MPa, except for a short transient the phase 
curve traces a limit cycle when the sand is not vibrated. The 
phase curve spirals towards a single fixed point, denoted by 
the filled green circle, when the material is vibrated. At the 
higher pressure of p = 50 MPa, except for short transients im¬ 
mediately following the switching on or off of fixed-intensity 
external vibrations, the phase curve traces a smaller limit cy¬ 
cle in the clockwise direction when the material is not subject 
to vibrations, and a bigger limit cycle when vibrations are 
switched on. 


pressure of p = 45 MPa, vibrations amplify the stick- 
slip amplitude before attenuating it, and steady sliding 
takes over at long enough times. Second, at the higher 
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pressure of p = 50 MPa, vibrations amplify the stick-slip 
amplitude and period. These observations are depicted 
in another manner in Fig. | 6 (b)] where we show the varia¬ 
tion of shear stress s with plastic strain rate 7 P* in a phase 
plot. The sudden onset or ceasing of vibrations results in 
a short transient in both cases, denoted by lighter shades. 
Subsequently the phase curve of s versus approaches 
and traces a limit cycle, or spirals towards a fixed point 
in the case of vibrated sand at p = 45 MPa. The limit cy¬ 
cle provides a straightforward illustration of the stick-slip 
amplitude. In other words, vibrations can either change 
the stick-slip amplitude or suppress stick-slip, depending 
on the other control parameters. 

Figure 0 shows that vibrations trigger the onset of slip, 
in the sense that vibrations clock-advance the first slip 
event at a reduced shear stress following the initial load¬ 
ing. This provides a possible explanation for dynamic 
earthquake triggering, the observation that seismic waves 
radiated by an earthquake can trigger catastrophic events 
elsewhere [ll|) [3 US ■ The triggering mechanism is as¬ 
sociated with weakening the threshold for slip events, 
rather than increasing the load. 


Vibrations advance stick-slip onset 



FIG. 7: (Color online) Shear stress s a s a function of accumu¬ 
lated shear strain 7, at applied shear rate 7 = 2.4 x 10^ s~^, 
with and without external vibrations at hxed intensity. The 
onset of slip occurs earlier in the presence of vibrations, as 
shown by the location of the first stress drop on the dashed 
curve for vibrated sand. This observation is reminiscent of 
of the triggering of an earthquake by seismic waves arriving 
from elsewhere. 

Finally, we show in Fig. | 8 ]our theoretical predictions 
regarding the pressure and shear rate regimes at which 
stick-slip operates. In addition to vibrations at inten¬ 
sity p = 5 X 10“^ as well as the unvibrated case studied 
earlier in the paper, we briefly examine here the effect 
of a higher vibration intensity p = 2 x 10 “^ on stick- 
slip behavior. (This corresponds to doubling the vibra¬ 
tion frequency, or quadrupling the amplitude. We set 
x{p = 2 X 10“^) = 0.16 to reflect the increased degree 
of compaction due to a higher vibration intensity.) We 
show that changing the vibration intensity p shifts the 
boundary of the parameter regime for which vibrations 



Imposed strain rate y (s *) 


FIG. 8: (Color online) Phase diagram showing the parameter 
space where stick-slip occurs, in the absence (p = 0), and 
at two different, fixed vibration intensities (p = 5 x 10“^^), 
2x 10~®. Stick-slip operates only at imposed shear rates above 
7 80 s“^ and pressures above p ~ 30 MPa. Vibrations 

suppress stick-slip for a small range of pressure at all shear 
rates for which stick-slip is possible (shaded in pink), but 
amplify the stick-slip amplitude for a wide range of shear rates 
and pressure (shaded in yellow). There is a narrow range of 
shear rates and pressure (shaded in magenta and green) for 
which stick-slip is possible in the absence of vibrations, and 
for only one, but not both, of the two vibration intensities 
p. The blue circles and red squares denote to the shear rate 
and pressure values corresponding to the stress-strain curves 


in Figs. [3] and 6(a) 


suppress stick-slip that would otherwise have been pos¬ 
sible. For our system of angular sand particles, there is 
a threshold normal stress Pmin ~ 30 MPa and threshold 
shear rate ymin ~ 80 s“^ below which only steady sliding 
is possible. Note from Eq. (EH) that the “effective spring 
stiffness” G/p is inversely proportional to the pressure, 
so the existence of a lower bound for the pressure corrob¬ 
orates with earlier results that indicate the necessity 
of a low spring stiffness for stick-slip instabilities. For 
each pressure value p above the threshold Pmin, stick-slip 
can occur only for a range of intermediate shear rates. 
This corroborates with Fig. [T] which suggests that insta¬ 
bilities are possible only for a range of imposed strain 
rates. 

Vibrations suppress stick-slip at the lower pressures at 
shear rates for which stick-slip would have been possi¬ 
ble. Of interest is the narrow range of shear rates and 
pressures, shaded in green in Fig. [51 for which stick-slip is 
possible only if the angular sand particles are unvibrated, 
or vibrated at the higher intensity p = 2x 10 “^, but that 
the lower-intensity vibration at p = 5 x 10“'^ would sup¬ 
press stick-slip behavior. The notion that some vibration 
amplitudes and frequencies can suppress stick-slip is cor¬ 
roborated by experimental observations and may be 
of practical importance to control engineering. At high 
enough shear rates, shear-induced dilatation described by 
xiq) dominates over the effect of interparticle friction ^ 
in Eq. (15.31) . As such, stick-slip becomes possible only 
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at very high pressures which, in addition to reducing the 
effective stiffness, also decreases the inertial time scale, 
and hence the dimensionless plastic strain rate q = 
as well as the amount of dilation described by 


VI. CONCLUDING REMARKS 

We have proposed a theoretical model for dense gran¬ 
ular flow. The model captures experimentally observed 
phenomena such as stick-slip instabilities. Additionally, 
it captures the clock advancement of slip due to exter¬ 
nal vibrations, accounting for the triggering of instabili¬ 
ties. The key ingredients of the model include: (a) noise 
produced by friction Eq. (I5.4|l . which increases the 
packing fraction when the granular material undergoes 
shear deformation; (b) a y-dependent yield stress pa¬ 
rameter ^ 0 , Eq. (14.13L which produces apparent rate¬ 
weakening behavior at intermediate strain rates; and (c) 
explicit account of work associated with external vibra¬ 
tions, Eq. dm. The underlying framework of nonequi¬ 
librium thermodynamics provides a natural way to in¬ 
corporate energy dissipation as well as the effects of in¬ 
terparticle friction and external acoustic vibrations. The 
compactivity, an effective-disorder temperature, plays a 
crucial role in the model. Its dynamics encapsulate how 
interparticle frictional dissipation and vibrations control 
the response of a granular medium to shear deformation. 

We found that noise generated by interparticle fric¬ 
tional dissipation constitutes an essential ingredient for 
stick-slip instabilities to occur, and found that this mech¬ 
anism is most prominent at intermediate shear rates 
(7 ~ 10^ s“^ for angular quartz sand) and high enough 
pressure. Additionally, we investigated the role of vi¬ 
brations in controlling stick-slip behavior; we found that 
vibrations trigger slip events, change the stick-slip ampli¬ 
tude and period over a range of pressure and shear rate 
values, and suppress stick-slip in other circumstances. 
Within our choice of material parameters, we have com¬ 
puted a stick-slip phase diagram that illustrates the con¬ 
ditions on the shear rate and normal stress that are nec¬ 
essary for stick-slip to occur. The fact that our model 
predicts experimentally observed phenomena, including 
auto-acoustic compaction [s^ and stick-slip instabili¬ 
ties, lends credence to the applicability of nonequilibrium 
thermodynamics to driven granular media. 

A physics-based model such as the one laid out by us 
in the paper offers a predictive description of phenomena 
that may occur at regimes not currently accessible by 
laboratory experiments. Our theoretical model offers a 
unified understanding of the underlying mechanisms for 
stick-slip failure, as well as amplification and suppression 
of catastrophic events. Based on physical principles as¬ 
sociated with grain-scale processes and nonequilibrium 
thermodynamics, these results shed light on practical 
problems in granular physics in earth science, industry, 
and beyond. For example, the notion that one can con¬ 
trol stick-slip behavior by external acoustic vibrations 


may have applications in the nondestructive evaluation of 
materials, the control of jamming and unjamming, dam¬ 
age, and wear in materials processing, as well as in other 
industrial applications. 

In the context of granular earth materials, an impor¬ 
tant example of phenomena observed in nature and in 
simulations, and captured in the present STZ model, is 
the dynamic triggering of catastrophic events [10, [ill, llO, 
[4^ . known to occur in the presence of a granular gouge 
layer. Specifically, we found that vibrations reduce the 
shear strength and clock-advance the first slip event upon 
loading. This corroborates with the observation that seis¬ 
mic waves emanating from one earthquake can trigger 
another elsewhere. Thus the STZ theory provides im¬ 
portant insight into the physics of fault failure in the 
earth in the presence of seismic waves. 

For much of this paper we have used one single vibra¬ 
tion intensity p corresponding to a single set of values of 
vibration amplitude and frequency. We see from Fig. [5] 
that a change in vibration intensity can control stick- 
slip behavior. To further characterize vibration-induced 
shear weakening and dynamic triggering, and to identify 
intervals of increased seismic hazard, it will be useful to 
investigate the effect of varying the noise strength p as 
well as the vibration amplitude and frequency. Ferdowsi 
et al. Eini have found in their simulations evidence 
of an amplitude threshold below which vibrations do not 
cause significant frictional weakening. It will be interest¬ 
ing to examine if the STZ theory can capture this obser¬ 
vation. This is beyond the scope of the present paper, 
but will be the subject of a future investigation. 

Our model results display evidence of preseismic slip, 
shown in Figs. |3] and [5] above. That is, the stress-strain 
curve shows some degree of rounding during the stick 
phase, corresponding to a gradual increase in the slip 
rate or plastic strain rate prior to the onset of slip 
events. Characterization of the stress state immediately 
prior to shear failure, as well as the degree to which pre¬ 
seismic slip hints at the arrival of large slip events, may 
have important implications on the evaluation of seismic 
hazards. 

In this paper we assumed for simplicity that the mate¬ 
rial is spatially homogeneous; in effect, we have consid¬ 
ered only what happens within the shear band. However, 
real granular materials are heterogeneous in all spatial 
directions, and the shear band thickness depends on a 
variety of factors such as boundary conditions, imposed 
shear rate, and aggregate shear modulus 0, jSO, IMl ■ 
It may be necessary to consider a spatially heterogeneous 
compactivity distribution in order to fully capture irr^- 
ular stick-slip dynamics seen in experiments such as [l|. 

We also assumed here for simplicity that grains do 
not break apart, in order to focus on the most salient 
physical ingredients associated with stick-slip. A simple 
mean-field physical model for grain fragmentation was 
discussed in earlier work [^ . There, we showed that 
grain fragmentation occurs readily above the crushing 
strength of particles, which is of the order of dozens of 
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MPa and determined partly by the grain size. Changes 
in particle size may have important implications on the 
frictional noise term the shear rate regime in which in¬ 
stabilities may occur, and the shear strength of the gran¬ 
ular medium. In we argued that grain fragmentation 
reduces the shear strength in a rate-strengthening gran¬ 
ular material. One could also easily envision the physical 
picture of small grains produced by fragmentation get¬ 
ting trapped between large grains, causing compaction. 
These aspects are beyond the scope of the present paper, 
but could have important implications on the stick-slip 
amplitude and frequency [2^. The effect of grain frag¬ 
mentation should be addressed in detail if one wants to 
examine the full range of possible physical behaviors in 
crushed rock particles. 

Future studies will account for the roles of inertial ef¬ 
fects, boundary conditions, the physics at the grain con¬ 
tacts, and grain fragmentation. These are in addition 
to addressing hysteretic phenomena and long-time relax¬ 
ation upon ceas ing acoustic vibrations or in unsheared 
granular media [l3, [HI [l^ , as well as an extensive char¬ 
acterization of stick-slip statistics such as dependence of 


the magnitude of stress drop and inter-event period on 
the vibration intensity p and other control parameters. 
While our results provide crucial input for refining the 
constitutive description of dense granular flow, we call 
for extensive experimental studies to validate our predic¬ 
tions and further constrain material parameters. 
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